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We study the dynamics of the scalar modes of linear perturbations around a flat, 
homogeneous and isotropic background in loop quantum cosmology. The equations 
of motion include quantum geometry effects and are expected to hold at all cur- 
vature scales so long as the wavelengths of the inhomogeneous modes of interest 
. remain larger than the Planck length. These equations are obtained by including 

holonomy corrections in an effective Hamiltonian and then using the standard vari- 
ational principle. We show that the effective scalar and diffeomorphism constraints 
■ are preserved by the dynamics. We also make some comments regarding potential 

00 ■ inverse triad corrections. 
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I. INTRODUCTION 



j> ■ The aim of loop quantum cosmology (LQC) is to use the methods and techniques of loop 

quantum gravity in order to study the dynamics of cosmological space-times when the space- 
! time curvature is of the order of the Planck scale and general relativity breaks down. For a 

| recent comprehensive review of LQC, see [1]. There has been great progress in LQC over the 

past few years, especially with regards to the simplest cosmological models, the homogeneous 
and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) space-times. The dynamics 
of the quantum theory have been studied extensively, using both analytical and numerical 
methods, and it has been found that the big bang singularity is resolved. In essence, quantum 
geometry effects introduce a strong repulsive force which causes a "bounce": there is a 
quantum gravity bridge between an earlier classical contracting universe and a later classical 
expanding universe [2, 3]. 

Now that the homogeneous sector of LQC is relatively well understood, the next step is 
to allow for small perturbations around a homogeneous and isotropic background and study 
how the presence of quantum geometry effects modify their dynamics when the space-time 
curvature is near the Planck scale. In this work we will focus on scalar perturbations; these 
modes are of particular interest in cosmology as they seed structure formation. With the 
results presented here, it will be possible to study the dynamics of scalar perturbations 
through the bounce for modes whose wavelengths remain larger than the Planck length. 

Since these results are of interest to researchers who are not experts in LQC, we have 
tried to ensure that this paper is as easy to read for nonspecialists as possible. If the reader 
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wishes to skip directly to the main result, the quantum-gravity-corrected equations for scalar 
perturbations around a flat homogeneous and isotropic background (presented in a notation 
familiar to cosmologists) are given in Sec. Ill B. 

In this paper, we will work with effective equations. Effective equations are obtained 
from a Hamiltonian constraint which has been modified in an appropriate way in order to 
incorporate quantum gravity effects. For homogeneous and isotropic models, effective equa- 
tions provide an excellent approximation to the full quantum dynamics of sharply peaked 
states, even at the bounce point when the quantum gravity effects are strongest [4], and we 
expect that this will continue to be true if small perturbations are allowed. 

There are two main types of corrections due to quantum gravity effects in LQC: holon- 
omy corrections and inverse triad corrections. Holonomy corrections arise as the connection 
is expressed in terms of holonomies of a minimal length and these corrections become im- 
portant when the space-time curvature nears the Planck scale. Inverse triad corrections 
are introduced in order to obtain well-defined operators corresponding to inverse powers of 
the area operator (this is necessary as is in the discrete spectrum of the area operator) 
and they play an important role when physical length scales become as small as the Planck 
length. We will only consider holonomy corrections here so that the equations hold for all 
curvature scales, but additional inverse triad corrections are necessary if any of the physical 
length scales of the space-time become comparable to the Planck length. 

Previous work studying cosmological perturbations in LQC at the effective level has fol- 
lowed a prescription, originally presented in Ref. [5], where the background and perturbation 
degrees of freedom are separated at the very beginning in the Hamiltonian constraint and 
the symplectic structure. This is certainly a possible approach to the problem, but it seems 
to be an unnatural way to try to incorporate nonperturbative quantum gravity corrections. 
For this reason, in this work we will separate the variables into the background and per- 
turbation parts only after the equations of motion have been determined. Once it is clear 
how to incorporate the holonomy corrections properly, it would of course be possible to go 
back and do it again using perturbed variables from the beginning, but the point is that the 
problem becomes less mysterious if one works with nonperturbed variables. 

Vector and tensor modes were the first to be studied using effective equations in LQC 
[6-8] . This was done using perturbed variables as in [5] and both holonomy and inverse triad 
corrections have been computed. Scalar modes have also been studied in [9, 10] where inverse 
triad corrections were included, and in [11, 12] where holonomy corrections were considered. 
Finally, there have also been some preliminary results about the phenomenology of quantum 
gravity effects that could potentially be detected in primordial gravitational wave signatures 
[13, 14] and in the cosmic microwave background [15, 16]. As the subtleties of the effective 
equations are better understood, phenomenological studies will become more robust and, 
we hope, will lead to falsifiable predictions. 

Scalar perturbations are the most interesting ones from an observational point of view 
and therefore it is important to understand how the equations of motion of the scalar 
perturbations are modified by quantum geometry effects. Since in homogeneous models, 
where the quantum dynamics are well understood [1-3], modifications due to holonomy 
corrections are significantly more important than those coming from inverse triad effects, 
we might expect this to continue to be true when small perturbations are allowed. This is 
why the study of holonomy corrections to the scalar modes of cosmological perturbations is 
important. 

While there has been some previous work studying holonomy corrections for scalar per- 
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turbations [11, 12], there remain many open problems. In Ref. [11], only large wavelength 
modes were considered and therefore all spatial derivative terms in the equations of motion 
were neglected. However, it is not clear that the wavelength of the modes observed today 
in the cosmic microwave background would necessarily have been large when the holonomy 
corrections are important (i.e., when the curvature nears the Planck scale). In some scenar- 
ios, such as inflation, this assumption clearly fails. Because of this, the results presented 
in [11] provide only partially corrected equations of motion for scalar perturbations. More 
recently, a further study has appeared where an effective Hamiltonian, including holonomy 
corrections, is studied [12]. However, the equations obtained from the effective Hamiltonian 
constraint are not consistent and one must add an ad hoc, non-unique modification to one 
of them in order to ensure that the constraints are preserved by the dynamics. It is not 
clear why the chosen modification should be preferred over any other. In this work, we will 
rectify the shortcomings of the previous studies of holonomy corrections to scalar pertur- 
bations by working with all wavelength modes and by introducing an effective Hamiltonian 
whose equations of motion automatically preserve the constraints. 

As mentioned above, we will work with nonperturbed variables in the effective Hamilto- 
nian constraint and symplectic structure. As we shall see, by doing this it will become clearer 
how to incorporate holonomy corrections from an LQC point of view. However, in order to 
work with nonperturbed variables in LQC, it is essential to have a diagonal metric. If the 
metric is not diagonal, a key mathematical simplification is lost: in homogeneous models, 
it is enough to only consider a small class of holonomies, called almost periodic functions. 
This is a great simplification which is not available in general, but can be generalized in a 
straighforward way to include space-times where the metric can be put in a diagonal form. 
Therefore, we will work in the longitudinal gauge where the metric is diagonal. This gauge 
choice is absolutely necessary if one only wishes to consider almost periodic holonomies: if 
one works in a different gauge (or without choosing a gauge), then one must work with a 
much larger class of holonomies and the problem becomes much more complicated. We will 
discuss this further in Sec. III. 

The outline of the paper is as follows: in Sec. II, we present a first order Hamiltonian 
formalism for cosmological perturbation theory which gives the standard equations for the 
scalar modes in general relativity. Then in Sec. Ill we modify the Hamiltonian in a suitable 
fashion in order to include holonomy corrections, this gives the quantum-gravity-corrected 
equations of motion for scalar mode perturbations. We also make some comments about 
inverse triad corrections in Sec. IV before ending with a discussion in Sec. V. 

In this paper, we will work in units where c = 1 but we will keep G and H explicit 
so that it will be possible to see whether a contribution is due to gravitational effects, 
quantum effects or both. We will define the Planck length as £p\ = y/Gh. Since we are 
only considering linear perturbations, all terms that are second order (or higher) in the 
perturbations will be dropped and thus all of the equations are understood to hold up to 
first order in perturbations. 

II. FIRST ORDER HAMILTONIAN FRAMEWORK 

In order to study small fluctuations around the flat FLRW cosmological background on 
a 3-torus, one allows small departures from homogeneity. A nice coordinate choice is the 
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longitudinal gauge in which case the metric is given by 1 

ds 2 = -(1 + 2^)dt 2 + a 2 (l - 2^)df 2 , (2.1) 

where a(t) is the scale factor and depends only on time while ip(x, t) encodes the fluctuations 
away from the mean scale factor and varies both with time and position. We have chosen 
the line element so that the volume of the 3-torus, with respect to the background metric 
given by ds 2 = da; 2 , is 1. Of course, one is free to make a different choice and one can check 
that this choice does not affect the results of this work. 

Let us justify the choice of the longitudinal gauge. This gauge choice is particularly useful 
as the resulting metric is diagonal and this considerably simplifies the situation as it will be 
possible to make a direct analogy with the homogeneous case without first perturbing the 
Hamiltonian constraint. In effect, a diagonal metric allows us to restrict our attention to 
holonomies that are almost periodic in the connection. It is precisely this simplification that 
has allowed the vast amount of progress in LQC over the past few years. This is also why it 
is difficult to include perturbations. Typically this means working with off-diagonal terms in 
the metric, which must either be treated in a perturbative manner, or one must consider a 
more general class of holonomies than almost periodic functions of the connection. However, 
as we shall see in Sec. Ill the longitudinal gauge allows us to avoid these two problems in 
the case of scalar perturbations. 



A. Elementary Variables 

In order to use a first order formalism, it is necessary to work with triads and co-triads 
instead of a metric. In the longitudinal gauge, the co-triads are given by 

w£ = a(l-#&i, (2.2) 

and the spatial metric is then given by q a } } = uj l a uu- The space-time indices a,b,c, . . . are 
raised and lowered by q ab while the internal indices i,j,k,... are raised and lowered by 
Sij = diag(l, 1, 1). Similarly, their inverse the triads are given by 

e? = 1(1 + ^)1?. (2.3) 

(X 

Note that since ^ < 1, we can drop all terms of the order of ip 2 or higher. The fiducial 
triads and co-triads used in the equations above are defined as 



In this work we consider the case of a massless scalar field. In this case (and for all other perfect fluids) 
there is no anisotropic stress, and therefore the variable <j) encoding the fluctuations around the lapse is 
necessarily equal to the variable ip which describes the deviations away from the scale factor. In some 
cases, this fails to hold in a quantum-corrected effective setting [10], but as we shall see the equations we 
obtain at the end of the following analysis are consistent and it does not appear to be necessary to allow 
for departures from <j> = ip here. 
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There exists a derivative operator D compatible with the triads and co-triads, 

D a e\ = d a e\ + T h ac el + ej?{e\ = 0, (2.5) 
where Y b ac is the usual Christoffel connection and Y l a is the spin-connection, 

K = -e* fc e$ (d [a u b]k + l -elu l a d [c u b] ^ . (2.6) 

The Christoffel symbols will not be used in this work, but it is necessary to calculate the 
spin-connection, 

Tl = -Tl = ed x ^, Tl = -Tl=ed y ^ rj = -T 2 X = ed z ^. (2.7) 

Here e = e 123 which can be ±1 and therefore e 2 = 1. 

Now, in order to study the perturbations in the (classical) loop gravity framework, it is 
necessary to use the Ashtekar connection and densitized triads as our basic variables. Since 
the metric is diagonal, we can parametrize the densitized triads by 

E^=p^e^ where p = a 2 (l-2^). (2.8) 

The parameter p is a function of position and time, but we will drop the arguments in order 
to simplify the notation, except where they are essential. 

The Ashtekar connection is given by A\ = T l a + where 7 is the Barbero-Immirzi 
parameter and K l a = K a i,e bl is related to the extrinsic curvature. It is easy to show that, 
like Ef, K l a is diagonal (with respect to the fiducial triads) and all of its entries are equal: 

K = {a- 2d^ - ai>)ul (2.9) 

A useful property of the Ashtekar connection in this context is that its diagonal terms solely 
come from K l a , while its off-diagonal terms solely come from T l a (of course, this is not the case 
in general). Since the densitized triads are diagonal, only the diagonal part of the Ashtekar 
connection will appear in the induced symplectic structure and therefore it is convenient to 
parametrize the Ashtekar connection by 

A x = cti - e(d z ip) t 2 + e(d y ip) r 3 , 

A y = e(d z ^j) n + cT 2 - e(d x il>) r 3 , (2.10) 
A z = -e(d y i/j) n + e(d x ^) r 2 + cr 3 , 

where A a = A l a Ti and the are a basis of the Lie algebra of SU(2) such that TiTj = 
\tij kr k — Following this definition, we find that the induced symplectic structure on 

our phase space gives the following nonzero Poisson bracket: 



{c(x),p(y}} = -^—5 3 (x-y). 



(2.11) 
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B. Massless Scalar Field 



Since we are primarily interested in the effects due to quantum gravity in the early 
universe, we will take the simplest matter field possible, a massless scalar field. The action 
for a massless scalar field is 



S = -\( vTiliT^) (2.12) 
and it is easy to show that the conjugate momentum to ip is given by 

n v = Ny/\ij\<p, (2.13) 

where the dot represents derivation with respect to time and N is the lapse function. The 
Poisson bracket is given by 

Mx),n v (y)} = 5 3 (x-y). (2.14) 

In this work, we will be using the Hamiltonian formulation of general relativity and 
therefore the stress energy tensor will not appear. Since a lot of the literature in the field 
of cosmological perturbation theory starts from the Einstein equations (and thus uses the 
stress energy tensor), the following dictionary can be useful in order to compare the results 
given here with those available in the literature. 

For a massless scalar field, the stress energy tensor is given by 

T"„ = g» x (d x p) (d u p) - ±8£ (Vpf • (2.15) 

For small perturbations around a flat FLRW background, we have tp — (p + 8tp where the 
bar denotes the homogeneous background and Stp is the inhomogeneous perturbation. One 
can immediately see that the second term in T fl u is negligible as it is second order in Stp. 

The degrees of freedom of a perfect fluid are given by the energy density p, the pressure P 
and the velocity four- vector . The energy density and pressure can be split into background 
quantities and perturbations, p = p + 8p, P = P + 8P, while the only scalar mode of the 
perturbation of the four-velocity at first order is given by 8u, which affects the spatial part 
of the four- velocity by 8u a = d a (8u). In terms of these variables, the stress energy tensor is 
given by (see, e.g., [17, 18]) 

T° = -p-8p, (2.16) 
T° a = (p + P)d a (8u), (2.17) 
T\ = SZ(P + SP), (2.18) 

and then by Eq. (2.15) the relations between the "standard" variables and the ones used in 
the Hamiltonian framework are the following: 

p = P = \p\ (2.19) 

8p = 8P = (p (8p) - (p 2 tp, (2.20) 
(p + P)8u = -06(p. (2.21) 
Another useful relation can be obtained by adding Eqs. (2.19) and (2.20), this gives 

P = P = (l-2V0y- (2-22) 
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C. The Hamiltonian Constraint and Dynamics 



Before writing the Hamiltonian constraint, it is necessary to express in terms of the 
conjugate variables (c,p). In order to do this, we will assume that f M ip = (i.e., the zero 
mode of the Fourier decomposition of ip is zero) and then it follows that 



p — p 



where p 



Jm 



(2.23) 



This requires the introduction of the nonlocal p, but it is necessary in order to express the 
perturbation ip in terms of the elementary variable p. As we shall see, the resulting nonlocal 
Hamiltonian system provides the usual local equations of motion for linear cosmological 
perturbations. 

The Hamiltonian constraint is given by 



C H = [ [NH + N a H a + N%], 
Jm 



(2.24) 



where, for a massless scalar field, the scalar constraint is 

1 



U 



16ttG7 2 



F ab «-(l + Y)V ab «) + 



7r* + E?E bi (d a <p)(d b <p) ) »0, (2.25) 



the diffeomorphism constraint is 

Ha 

and the Gauss constraint is 



Fab + K^daip « 0, 



ft = V a E? = d a E? + e^A^El » 0. 
The F ah k appearing above is field strength of the Ashtekar connection, 

F ab k = 2d [a A k b] + e i3 k A\A{. 



(2.26) 



(2.27) 



(2.28) 



Similarly, Q ab k is the field strength of the spin-connection. 

The ~ indicates that T-L,T-L a and Qi are constraints and must vanish for physical solu- 
tions. A solution to these three constraints at an initial time can then be evolved to later 
times by taking Poisson brackets with the Hamiltonian constraint Ch- An important point 
is that any solution that initially satisfies the scalar, diffeomorphism and Gauss constraints 
will continue to do so under the evolution determined by Ch- 

For linear perturbations around the flat FLRW space-time in the longitudinal gauge, we 
find that the Gauss constraint is automatically satisfied, while the scalar constraint becomes 



n 



\p\ f 3c 



8nG 



r 



2V 2 (^ 



2p 



and the diffeomorphism constraint is given by 



"Ha 



d a c + cd, 



i ( PzJp\ 

a \ 2p J 



2p 3 / 2 q 2 \ * 



+ TT^CV/? ~ 0. 



0, (2.29) 



(2.30) 
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Finally, since N — 1 + — 1 + ^-f and iV a = in the longitudinal gauge, the Hamiltonian 
constraint is given by 

Ch= Li l+p: ^) n ' (2,3l) 

Although one might be tempted to drop terms in % like (V(p) 2 which are second order in 
the perturbations, they contribute first order terms in the evolution equations (in this case 
for -ftp) after the derivative is integrated by parts. 

This Hamiltonian formalism is of course well-defined for all c(x),p(x), (p(x) and 7i v (x), but 
the only case we are interested in is a homogeneous background with small perturbations in 
which case the resulting system is equivalent to linear cosmological perturbations in general 
relativity. Therefore, we will only consider equations for linear perturbations around a 
homogeneous background in what follows. 

It is now possible to obtain the standard results for linear perturbations in cosmology from 
the two constraints and the Hamiltonian. Using the relations (but only after the equations 
of motion have been obtained by the variational principle) 

p = a 2 (l - 2*0), p = a 2 , c = c + 5c, <p = <p + 5<p, — + Sn^,, (2.32) 

where as usual the "barred" quantities correspond to the unperturbed background quantities, 
we find that the scalar constraint H — implies that 

ac 2 , 2<2 9 , AttG vt 2 

— O-0--S-VV+— ^j = 0, (2.33) 

where we have only kept terms linear in the perturbations. The diffeomorphism constraint 
Ha = gives 

— \d a (5c) + cd a ijj] + AnG^d a {5ip) = 0. (2.34) 
7 L J 

The dynamics are obtained from the Hamiltonian constraint. Starting with the equations 
of motion for the matter degrees of freedom since they are simpler, 

5C h _ . n^2, 



^ = {k v ,C h } = = VqaV 2 (5if). (2.36) 

Note that the first relation is equivalent to Eq. (2.13), as one should expect and it also 
implies, together with Eq. (2.22), that 

P = 7^L, (2-37) 

which will be a useful relation later. 

The equation of motion for p also has a simple form, 

V = {P, C H } = = i (1 + ^ )VWh 

3 be 7 
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which, solving for c, gives 

aij?j , (2.39) 

which is exactly what is given in Eq. (2.9) (recall that the diagonal portion of A l a is given 
by 7 fQ as the diagonal part of T l a is zero). 

Finally the equation of motion for c is the most complicated as one must vary Ch with 
respect to p. There are two contributions that must be integrated by parts (since the 
topology is T 3 there are no boundary terms) and there are also terms which depend on 
p which could potentially contribute. Since 8p/5p(x) = 1 — note that there is no delta 
function — all of the terms that are obtained by varying p with respect to p(x) are inside an 
integral over Ai and would contribute a nonlocal term to the dynamics; this would disagree 
with the well known general relativity results. However, all of the nonlocal terms vanish and 
one is left with a local equation of motion. In order to see this, we start with 



c(x) = {c(x),Ch} = ^- I d 3 y 

6 Jm 



(2.40) 



where the first term can be dropped as it is multiplied by the constraint H which vanishes 
on-shell. A straightforward calculation then gives 

= _ [l+^(x)]c(x) 2 _ ^ m 2vr 7 G[l + 4>(x)]n,(xy 



2 7v /p(£) 3a rv p{xfl 2 

The fourth term can dropped as it is a total divergence, while the fifth term requires a little 
more work. From Eq. (2.23), it follows that 

fy(v) = K$ _ 63 ^ Zll (2 42) 

Sp(x) 2p 2 2p ' 1 ' ' 

where the delta function only appears in one of the two terms. Using this relation, integrating 
by parts and dropping terms that are second order in ip, the fifth term in Eq. (2.41) becomes 



^ J/y[vm\ v (||) = - 1 / dfW(fi). (2-43) 



The second term vanishes as it is a total divergence and one is left with the local relation 

1 7T 2 

6= - — (1 + 2^)c 2 - 27rG 7 a^. (2.44) 
2 7 a 

In order to compare our results to the standard results in the literature, one can remove 
all instances of and c via Eqs. (2.35) and (2.39) and then use the definitions of bp and 
8u given in Eqs. (2.20) and (2.21). After doing this, the scalar constraint gives 

d 2 8ttG 8ttG d 2 d • 2 2 . . 

- 2 = — p + — Sp + 2-* + 2-J - - 2 VV, (2.45) 
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and the diffeomorphism constraint % a = becomes 

d a (dip + aij?j + 8nGapd a (Su) = 0. (2.46) 
Equation (2.36) multiplied by tp gives 

p + 6-p + 5p + 6-5p - 6pij + ^pV 2 (5u) = 0, (2.47) 
a a a z 

while an additional relation between the perturbations in the massless scalar field can be 
obtained by taking the time derivative of Eq. (2.21) which gives 

d t (2p5u) = -bp - 2pip - 6-pSu, (2.48) 

(X 

where the relation (p = — 3a<p/a — obtained from Eqs. (2.35) and (2.36) and the relation 
p = 2aa — was used. Finally, using Eq. (2.39), Eq. (2.44) gives 

... 2 • 2 

- + At = -AitGp - AitGdp + 2-i\) + + 4-i> + (2.49) 
a 2a 2 a a 2 a 

These five equations provide the standard results for linear perturbations around the flat 
FLRW background (compare with, e.g., [17, 18]) when the matter content is a massless 
scalar field; that is, when there is no anisotropic stress and the equation of state is P = p. 
This shows that the Hamiltonian (2.31) with the constraints (2.29) and (2.30) gives the 
correct dynamics for scalar perturbations in cosmology. 



III. HOLONOMY CORRECTIONS 

In this section, we will appropriately modify the scalar constraint in order to incorporate 
quantum gravity effects that come from using holonomies of the connection — rather than 
the connection itself — in the quantum theory. On the other hand, the diffeomorphism and 
Gauss constraints are not modified: this is because they encode symmetries which we expect 
to remain present at all scales, including the Planck scale. 

The idea is to model the most important quantum gravity effects coming from a quan- 
tum theory in a simpler (often called "quantum-corrected") setting which is treated in a 
classical manner. One example of a quantum theory that could generate these quantum- 
corrected equations is a form of lattice LQC where each cell is described by a homogeneous 
and isotropic geometry but where the geometries vary from one cell to another so that in- 
homogeneities are present at scales larger than the cell size [19, 20]. In this setting, the 
Hamiltonian would be separated into an ultralocal term (the "homogenous" part) and all 
of the terms containing derivatives would be treated as interactions between different cells. 
This would of course constitute an approximation but would be valid so long as the volumes 
of the individual cells remain larger than the Planck volume and that there are enough cells 
in the lattice in order to properly describe the inhomogeneities up to a prescribed maximum 
wave number. 

In the longitudinal gauge, the homogeneous metric in each cell is isotropic and thus the 
symplectic structure in each cell would be that of an FLRW space-time. This indicates 
that one might be able to quantize each cell following the methods used in homogeneous 
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models [1-3] and then turn on interactions in the Hamiltonian. It is worth pointing out 
that the diagonal form of the metric in the longitudinal gauge will greatly simplify the 
quantization process. This is mainly because, as one can easily check, when the metric 
is nondiagonal the field strength obtained by taking a holonomy around a square loop of 
nonzero area is generically not almost-periodic in the connection (which is proportional to 
the extrinsic curvature in the spatially flat homogeneous case where the spin-connection 
vanishes). However, in this lattice LQC quantization, due to the diagonal form of the metric 
it is possible to restrict our attention to almost periodic functions of the connection, just as 
in homogeneous space-times. 

Once the quantum theory has been constructed, the effective theory can be obtained by 
taking a sharply peaked wave function and studying the evolution of the expectation values 
of certain observables under the action of the Hamiltonian operator. Some observables such 
as p can be obtained directly from expectation values of the fundamental operators in the 
quantum theory while others — such as c which one could associate to the expectation value 
of an operator corresponding to a short holonomy (minus the identity) divided by the length 
of the holonomy — would require a little more care. In any case, we will not construct the 
quantum theory here but leave it for future research. Instead, we will build the effective 
Hamiltonian by applying the heuristic methods following from (and well supported by) 
homogeneous LQC. 

In the first part of this section, we will implement holonomy corrections in the system 
studied in the previous section, that of a massless scalar field. Then, in the second part we 
will present a conjectured generalization of these results for all perfect fluids. 

A. Effective Equations for a Massless Scalar Field 

One of the key parts in the process of the LQC quantization of the Hamiltonian constraint 
operator is that the field strength F ab k of the Ashtekar connection must be expressed in terms 
of a holonomy of the connection A % a around the smallest possible loop. It is posited that 
the smallest loop has an area given by A£ P1 = 4y / 37r7£p 1 , the smallest eigenvalue of the area 
operator in loop quantum gravity. 

For homogeneous and isotropic cosmologies, it has been found that the full LQC dynamics 
of states which, at some initial time are sharply peaked around a classical solution, remain 
sharply peaked as they are evolved by the action of the Hamiltonian constraint operator [2]. 
In addition, the dynamics of the wave packet are extremely well approximated by a modified 
classical Hamiltonian constraint where one includes holonomy corrections by determining 
the field strength by taking a holonomy around the minimal loop as prescribed above [4] . In 
this sense, it seems that the most important quantum gravity modifications to the Einstein 
equations in this very simple setting come from the holonomy corrections. These corrections 
become important when the space-time curvature approaches the Planck scale, but are 
completely negligible when the matter energy density is less than O.Olppi. 

Unfortunately, a priori it is not entirely clear how to incorporate holonomy corrections 
for inhomogeneous space-times. In previous works [6-8, 12], generic holonomy correction 
terms were added to the effective equations by hand and then the form of the correction 
terms was restricted by the condition that the constraint algebra must continue to close. 
However, so far it has proven to be impossible to complete this programme for holonomy 
corrections to scalar perturbations. 

Here, due to the choice of both the gauge and of the variables, the situation is luckily quite 
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simple. In the classical scalar constraint, we find that the term with the difference between 
the field strengths of the Ashtekar connection and the spin-connection can be simplified, 



E«E b ^ k {F ab k - (1 + 7 2 Rb1 = E^E b /h[Ft )k ~ ^J] 



(3.1) 



where we have defined 



^(iso)fc 
ab 



(3.2) 



where of course c depends on position. Although this "field" strength is also inhomogeneous, 
it nonetheless has a very similar form to the field strength in homogeneous settings. Note 
that this provides a natural separation of the F ab k and Q ab k terms into an ultralocal term, 
F^°^ k , and an interaction term, Vt ab k . Therefore, following the arguments introduced at 
the beginning of this section, we can introduce holonomy corrections to the ultralocal term 
^ab°^ k as m homogeneous models while Q ab k is treated as an interaction term which only 
depends on p and therefore will not be modified. Thus, in order to incorporate holonomy 
corrections in the scalar constraint, we shall mimic the procedure in the homogeneous case 
and replace 2 



^(iso)fc 
ab 



c 2 e k ij u i a uji -> 



s_mjxc k 

— 2 6 ijUcU b , 

fi 



(3.3) 



where 



\P\ ' 



just as in the homogeneous case [1, 2], except that p 
Thus, the scalar constraint becomes 



(3.4) 

a 2 (l — 2ip) is no longer homogeneous. 



H 



cff 



-1 f3\p\ 3 / 2 2 



8nG V7 2 Af 



+ 



Ti Z 



2\p\ 3 / 2 q 



+ 



cos 



2jxc (vv?y 



(3.5) 



where the cos2/2c in the last term is added in order to ensure that the effective scalar 
constraint is preserved by the dynamics, i.e., that "H cff ~ 0. 

On the other hand, the diffeomorphism constraint is not modified: 



n 



cff 



Vqp 



d a c + cd a [^ 



TTpdaif PS 0. 



(3.6) 



2 The Belinskii, Khalatnikov, Lifshitz (BKL) conjecture — which suggests that spatial derivatives can be 
ignored when the spatial curvature becomes very large [21, 22] — also supports this way of introducing 
the holonomy corrections. Since holonomy corrections are only important when the curvature approaches 
the Planck regime, if the BKL conjecture holds near the Planck scale then spatial derivatives and thus 
£l a b k can be ignored when one incorporates holonomy corrections. Then only F^°^ k remains and the 
holonomy corrections are introduced as in Eq. (3.3). Note that although the BKL conjecture says nothing 
about isotropy, F^°^ k is locally isotropic due to the locally isotropic form of the metric in the longitudinal 
gauge. 
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The Gauss constraint is not modified either and continues to hold due to the choice of the 
longitudinal gauge and therefore we do not need to consider it any further. As an aside, we 
point out that it is not immediately clear how to incorporate an unchanged diffeomorphism 
constraint in lattice LQC, we leave this problem for future work. 

As before, the Hamiltonian constraint provides the dynamics and is simply given by 

C* = J (l + fc?) H cS . (3.7) 

The scalar and diffeomorphism constraints must be satisfied by the initial data and thus 

\p\ 3/2 9 2a „ 9 , 4nG 9 , . 

■ sin 2 Jic - — V 2 ^ + ^ £ = 0, (3.8) 



7 2 A£ P1 r 3 Y 3\p\ 3 / 2 q * 

V^^^ [9a(5c) + c9 a ^] + n v d a (5 V ) = 0. (3.9) 
Finally, the equations of motion generated by Cff are the following: 

* = < 1 + *>^' (310) 

7r v = A/^acos2/ic V 2 (/?, (3-11) 

p = = sin lie cos //c, (3-12) 

7V A£pi 

3a . 2 _ c(l + -0) . _ _ 27rG7 2 

c = ~~ - — r W sm ^ c ^ — sin //c cos //c g — (1 + bipjn^. (3.13) 



a" 



Once again all nonlocal contributions due to a variation of p vanish, as one should expect. 
Note that while classically c is proportional to p, this is not the case in the effective theory 
and the relation between them, as seen in Eq. (3.12), is considerably more complicated. 
Finally, the relationship between ti^ and p given in Eq. (2.37) continues to hold. 

It is interesting to solve for the Hubble rate by squaring Eq. (3.12) and using Eq. (3.8), 
this gives the modified Friedmann equation (including perturbations), 

where the critical energy density is p c = 3/(87rG'7 2 A£p 1 ). Clearly, the classical result is 
obtained so long as the local matter energy density remains well below the critical energy 
density. It is also clear that when the matter energy density approaches the Planck scale, 
the Hubble rate diminishes due to the repulsive nature of the quantum gravity corrections, 
reaches zero close to the critical density (depending on the local strength of the fluctuations 
in ip) at which point there is a bounce where the Hubble rate begins to grow again and 
as soon as the matter energy density drops below the Planck regime, contributions due to 
quantum gravity are completely negligible. 

Finally, the form of the equations of motion given above is a little unwieldy and for some 
applications it will be useful to expand these equations and separate the background terms 
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from the perturbed quantities. For the scalar field, this is easy, 

(M = ^A ("SI 

n v = 0, (3.17) 

(6n v ) = A/gacos2/icV 2 5y?; (3.18) 

but it becomes a little more involved for the geometric degrees of freedom. For p, since 
p = a 2 (l — 2*0), we find that 

It is also necessary to expand the trigonometric functions of fic. This gives, for example, 

sin uc = sin 1 c0 + be cos . (3.20) 

a a a 

From these equations, it follows that 
. . a 2 . 2y^f P1 c 

p = 2aa = — t= — sin , (3-21) 



pi 



a 



• -0 . 2>/A£pic ap + Sc 2VA£ Pl c 

tp — = — sin cos ; (3.22) 

27VA£pi a 7a a 

-3a . 2 ^pic , c 2y/A£ P1 c 2nG-y _ 2 , 

sin 1 p= — sin — 7r , (3.23) 



27A£ P1 a 2 7 \/A£ 



pi 



a a° 



c(cip + 5c) 2y/A£ Pl c c^j + 5c . 2y/A£ Pl c UttG^ _ 2 47rG 7 _ r 

(5c) = cos ? =— sin vr — n v 6n v . 

7a a 7v / A£ P1 a 5 ^ a 5 

(3.24) 

The trigonometric identities 2 sin cos = sin 20 and cos 2 — sin 2 9 = cos 20 have been used 
in order to shorten the expressions above. 

The dynamics, both of the background and of the perturbations around it, are given by 
the equations above and one can check that, on-shell, H cS = Uf = 0, just as one would 
expect. Note that the cos2/2c term that was added to the effective scalar constraint in Eq. 
(3.5) is essential for "H cff = to hold. This result shows that the equations of motion and 
the constraints are consistent with each other. 



B. Conjectured Generalization to Other Perfect Fluids 

In this part, we provide a conjectured extension of our results to all other matter fields 
that can be treated as perfect fluids. It would be possible to use a Hamiltonian framework 
for perfect fluids [23], but we will show that this is not necessary: the generalization is 
straightforward enough that it can be done by hand. There is a possibility that a strong 
quantum backreaction from the matter field could spoil the validity of the effective equations, 
but we will not consider this potential effect here. 
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It is worth pointing out that it would be significantly more difficult to generalize these 
equations in order to include matter fields that allow anisotropic stress since at the very 
beginning it was assumed that the perturbation in the lapse, 0, was equal to the perturbation 
in the scale factor tp. If this assumption is not made, then it is not immediately clear how 
to construct a Hamiltonian formulation with the symplectic structure given in Eq. (2.11), 
without expanding it in terms of the background and perturbations which would ruin the 
simple substitution used in Eq. (3.3). 

Nonetheless, this generalization to all perfect fluids will be very useful, especially since 
many matter fields behave like radiation at very high energies, which are precisely the 
conditions to be expected when the curvature is near the Planck scale. 

For generic perfect fluids, the effective scalar constraint becomes 



P ~ 



8ttG 



-PI 



3a 2 



o, 



(3.25) 



which, as sin 2 9 < 1, shows that the matter energy density is bounded above by the critical 
density 

Pc = (3 ' 26) 

(modulo some small corrections from V 2 ip), which is consistent with Eq. (3.14). 
As before, the diffeomorphism constraint is unchanged and therefore 



1 r 

7 



cd a i) + d a (6c) - ArrGa(p + P)d a (5u) = 0. 



(3.27) 



Note however that since c is not proportional to p when holonomy corrections are present, 
the diffeomorphism constraint cannot be simplified as much as it usually is when there are 
no holonomy corrections. 

The equations of motion for the matter degrees of freedom are easy to generalize. Using 
Eqs. (2.37) and (3.11), along with the scalar constraint in order to solve for the cos2/ic term, 
we find that 



p + 3-(p + P) + 5p + 3-(5p + 5P) - 3(p + P)j) + — (p + P)V 2 (8u) 1 



2p 



0, (3.28) 



where the critical density p c appears once more. The equation for Su, 



(p + P)5u =-5P-(p + P)4>-3-( y p + P)5u, 



(3.29) 



continues to hold. 

Moving on to the geometrical degrees of freedom, since there are no matter terms in the 
equation for p it is trivially generalized, 



. 2p(l + ^) . _ 

p = -= sin pc cos uc, 



(3.30) 



while the generalization of Eq. (3.13) is a little trickier. A careful analysis of the equation 
of motion for c at the classical level shows that the 7T 2 , appearing there is related to the 
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pressure of the massless scalar field rather than its energy density. Therefore, it follows that 
for an arbitrary perfect fluid, 

3a • 2 - . C ( X + V0 • - - a n v> /q qi \ 

sin jic H = — sm jiccosjic — AnG^aF. (3.31) 



2 7 A4, ^ 7v /A£ 



pi 



Depending on the situation, it may be useful to expand these equations as was done in the 
previous subsection in Eqs. (3.15) — (3.24); this can be done by following the same procedure 
that was used there. Also, one can check that for these effective equations, generalized 
in order to include all perfect fluids, the scalar and diffeomorphism constraints are again 
preserved by the dynamics. 

We can obtain the modified Friedmann equation by squaring Eq. (3.30) and using Eq. 
(3.25), 

d 2 8ttG / p\ h 2 , h ■ 2 - , / 2p\ , n nn . 

7 2 = — p ( 1 " J ) + 2 ^ + ^ ~ 3^ V V ( 1 - 7 ) ' (3 ' 32) 

LI O \ Uq /Lb LI OLl \ Uq / 

where it is understood that p = p+Sp and the critical density p c appears again. Similarly, the 
modified Raychaudhuri equation can be obtained by taking the time derivative of Eq. (3.30) 
and then using Eqs. (3.25), (3.30) and (3.31) in order to remove all of the trigonometric 
functions. This gives 



a a 2 'a . h 2 , a ■ 

- = — + 2-ij-2—ij + -ij + ij- 

a a A a a z a 



4ttG(p + P)-—V 2 ^ 



!-- + ^7^V^) (3.33) 
p c 2nGp c a 2 



These equations can be compared with the usual general relativity results obtained in Sec. 
IIC, the classical results are obtained in the limit p c — > oo. 

Any initial data satisfying Eqs. (3.25) and (3.27) can then be evolved by using Eqs. 
(3.28), (3.29), (3.32) and (3.33); these last four equations are written in a form similar to the 
standard cosmological notation. Although the constraints are written in terms of different 
variables, one can use the standard classical constraints given in Eqs. (2.45) and (2.46) in 
order to constrain the initial data so long as the energy density p is far from the critical 
density p c at the time the initial conditions are set. If the initial conditions are chosen at a 
time where p is within two orders of magnitude of p c or less, then the holonomy-corrected 
constraints (3.25) and (3.27) must be used. 



IV. SOME COMMENTS ON INVERSE TRIAD CORRECTIONS 

In this section, we will make some comments regarding the type of inverse triad corrections 
that might be expected in the effective equations for cosmological perturbations. This section 
stands apart from the remainder of the paper and can be skipped on a first reading. 

There are two main inputs from loop quantum gravity that are used in order to con- 
struct the Hamiltonian constraint operator in LQC. First, the field strength of the Ashtekar 
connection must be expressed in terms of holonomies and second, inverse triad operators 
must be introduced carefully as the operator p includes zero in its discrete spectrum. In the 
previous section, we showed how it is possible to treat holonomy corrections in an effective 
action; in this section we shall comment on some of the properties that we expect inverse 
triad corrections to have. This has already been studied in LQC (see, e.g., the appendix 
of [24] where similar considerations arise in a slightly different setting) but there are some 
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important new aspects which appear in the context of lattice LQC which provide further 
insights into this problem. 

As is well known, there are many different ways to build inverse triad operators in LQC, 
all of which have the correct semiclassical limit but behave differently at the Planck scale. 
Because of this ambiguity, it is not clear which inverse triad operator to choose. Although 
there are some heuristic arguments indicating that inverse triad corrections might be com- 
parable to holonomy corrections [15], all of the inverse triad operators defined in LQC so 
far vanish in noncompact space-times. This indicates that our current understanding of 
inverse triad operators is flawed as these operators should be local, rather than global, op- 
erators. However, the inverse triad operators that naturally arise in the setting of lattice 
LQC have different properties and may explain what the correct way to implement inverse 
triad operators in a cosmological setting is. 

In order to gain an idea of the form inverse triad corrections should take, we will consider 
a relatively simple lattice LQC model [19, 20]. In this setting, we discretize the 3-torus into 
A 3 cells which are each taken to be homogeneous, but as the gravitational and matter fields 
can vary from one cell to the next, inhomogeneities are present at large scales. The number 
of cells, A 3 , gives the largest wave number perturbations may have in this discretization, 

A 

k=-. (4.1) 

Since the gravitational field is taken to be homogeneous in each cell, the induced symplectic 
structure from the full theory changes and the Poisson brackets become 

87^7 A 3 



{c(n),p(m)} 



(4.2) 
(4.3) 



{<p(n), 7i>(m)} = N 3 5 H)Tfl , 

where n,m are vectors which label the cells in the three-dimensional lattice and is a 
Kronecker delta. Note the explicit dependence of A 3 in the symplectic structure, this will 
reappear in the inverse triad operators. 

Following the full theory, inverse triad operators are constructed by starting with a Pois- 
son bracket (which is equal to 1/p when it is evaluated) between two functions on the phase 
space that can easily be promoted to operators and then the inverse triad operator is given 
by the commutator between the two operators, divided by ih [25]. For lattice LQC, a simple 
choice for the inverse triad operator corresponding to 1/p acting on a particular cell is given 
by 3 



p(n) 



\p(n)) = 



47r 7 v / A4 1 A 3 



|p(n)| 3/2 + 2 7 r 7 v / A4 1 A 3 



1/3 



\p{n)\ 



3/2 _ 



27T7V A^jA 3 



1/3 



\p(n)). 



(4.4) 



3 This can be derived by noting that p- 1 = -(3i/4TrGj\^e P iN 3 )e-^ c / 2 {e^ c , y^\p\}e- 1 ^ 2 (where we 
have suppressed the label n in order to simplify the notation) and then promoting the right side to be 
an operator by replacing {•, •} — > [•, -]/ih. The action of this operator is most easily seen by changing 
variables to v oc p 3 ' 2 so that e l ^ c acts as a simple shift operator on states \u). Then, changing variables 
back to p one obtains the result given in Eq. (4.4). See [1, 2] for further information on inverse triad 
operators in the p, scheme in homogeneous and isotropic LQC. 
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Clearly this operator is always well-defined, even on the state \p(n) = 0), which it annihilates. 
Also, the eigenvalue approximates 1 jp for large p. 

In order to incorporate inverse triad corrections for the operator chosen above, one would 
simply replace all occurances of 1/p in the scalar constraint with the right hand side of Eq. 
(4.4). Since we know that the singularity is avoided due to the holonomy corrections, it 
follows that p 7^ and then inverse triad corrections can be written as 

- ^- x F(p,N), (4.5) 
p P 

where the corrections are encoded in the function F(p, N), 

F(P,N) = [\l + /(p, 7V)|V3 - |1 - f(p,N)\^} , (4.6) 

and f(p, N), in turn, is given by 

^^m. (4 . 7) 

Since the inverse triad correction depends on the number of cells N 3 , this type of cor- 
rection is sometimes viewed as being unphysical. However, since the physical wavelength of 
the Fourier mode in a compact, (approximately) homogeneous, space of volume p 3 ^ 2 divided 
into iV 3 cells is A p h y = 2^/p/N [see Eq. (4.1)], it follows that F and / should not be viewed 
as functions of p and N but rather of the physical wavelength of the Fourier mode studied, 
Aphy. Then, 

^(A phy ) = [|1 + /(A phy )| 1/3 - |1 - /(A phy )| 1/3 ] , (4.8) 

7(A phy ) = (4.9) 
Vy 

The inverse triad correction to a given Fourier mode depends on the ratio of the wavelength 
of that mode to the Planck length. For modes with a wavelength considerably larger than 
£pi, inverse triad corrections are completely negligible, even if the space-time curvature is 
large. 

The dependence of the inverse triad correction on the wavelength of each mode is com- 
pletely natural. In homogeneous models in LQC, inverse volume effects are only relevant in 
compact space-times and then the strength of the correction depends on the ratio of the vol- 
ume of the space-time to the Planck volume. This indicates that for inverse triad effects to be 
present, there must be length scales in the space-time which can be compared to the Planck 
length. In homogeneous space-times there is only the size of the entire space-time which 
provides a length scale 4 , but when perturbations are present there are additional length 
scales provided by the wavelengths of the Fourier modes of the perturbations. Therefore, 
this property of the inverse triad corrections should not be surprising. 

As a final note, we point out that much of the previous work studying inverse triad 
corrections in the cosmology perturbation equations for scalar, vector and tensor modes do 
not allow for the inverse triad corrections to depend on the wavelength of the mode, instead 
the corrections typically only depend on the volume of the background homogeneous space- 
time. We suggest that these results should be generalized. 



This is why inverse volume corrections in homogeneous, noncompact space-times vanish in LQC. 
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V. DISCUSSION 

We have shown how it is possible to incorporate holonomy-type corrections into the 
equations of motion for perturbations around a flat FLRW background. This was first done 
in the case of a massless scalar field and then we presented a conjectured generalization for 
all perfect fluid matter fields. 

Since holonomy corrections are nonperturbative, we decided to use elementary variables 
that are nonperturbative. This is to say that they are not separated into background and 
perturbation terms, rather it is after the equations of motion are obtained that this split is 
performed. It was then necessary to work in the longitudinal gauge so that we could restrict 
our attention to holonomies that are almost periodic in the connection. These two choices 
greatly simplified the problem and the correct way to implement holonomy corrections in 
the effective Hamiltonian constraint for this particular gauge became more obvious. 

It would be nice to extend this analysis in order to include tensor and vector modes, but 
this may not be as trivial as one would hope. The problem is that a lot of the simplifications 
which made this problem tractable at a nonperturbative level — such as a diagonal metric 
which allowed us to work with a relatively small phase space and restrict our attention to 
holonomies that are almost periodic functions of the connection — are not possible once 
tensor and vector perturbations are allowed. Therefore it seems that a combination of 
the techniques used here and perturbative methods used in [6-10, 12] will be necessary in 
order to be able to treat scalar, vector and tensor perturbations all together. Despite these 
remaining open issues, a lot can be learnt from this treatment of scalar perturbations. 

First, the diffeomorphism and Gauss constraints are unchanged and therefore the symme- 
tries they encode hold exactly at all scales, including the Planck scale. On the other hand, 
the scalar constraint and the dynamics of the gravitational field are modified by the holon- 
omy corrections. As we saw in Sec. Ill, the holonomy corrections are introduced in precisely 
the same way as in the homogeneous models and the resulting constraints are shown to be 
preserved by the dynamics. As usual in LQC, we found that the holonomy corrections are 
parametrized by the critical energy density p c ; these corrections vanish in limit of p c — > oo 
and then the standard cosmological perturbation equations are recovered. 

Although it is not yet clear how to include inverse triad operators — which are essential 
if one wants to construct the quantum theory — it is possible to learn some lessons from 
a discretization of LQC. The main point is that the inverse triad corrections should vary 
from one Fourier mode of the perturbation to the next. To be more precise, the inverse 
triad corrections, for a given Fourier mode of the perturbation, should explicitly depend 
upon the ratio of the wavelength of the mode to the Planck length. Thus, the shorter the 
wavelength, the more important the inverse triad correction will be. This is natural as 
the only physical length scale provided by perturbations is the wavelength of the modes. 
Because of this, inverse triad corrections should be expected to be completely negligible so 
long as the wavelengths of the perturbative modes of interest stay an order of magnitude or 
two away from the Planck scale. However, if the wavelengths of interest approach the Planck 
scale, then inverse triad corrections cannot be ignored. For example, these corrections could 
be important in inflationary scenarios, but not in ekpyrotic or matter bounce models. 

This shows that holonomy corrections and inverse triad corrections have a completely 
different nature: holonomy corrections become important when the space-time curvature 
approaches the Planck scale while inverse triad operators become important when physical 
length scales are comparable to the Planck length. Thus, depending upon the situation, one 
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of the corrections may be important while the other is negligible. In particular, we expect 
the equations derived in this paper hold at all curvature scales so long as the wavelengths 
of the scalar modes of interest remain large compared to £p\. 
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